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ABSTRACT 

We present the largest sample to-date of intermediate-resolution blue-to-red optical spectra of B-type su- 
pergiants in M 3 1 and undertake the first survey of diffuse interstellar bands (DIBs) in this galaxy. Spectral 
classifications, radial velocities and interstellar reddenings are presented for 34 stars in three regions of M 3 1 . 
Based on a subset of these stars with foreground-corrected reddening E^}y > 0.05, the strengths of the M31 
DIBs are analysed with respect to the amount of dust, the ultraviolet radiation field strength and the PAH emis- 
sion flux. Radial velocities and equivalent widths are given for the /i5780 and /16283 DIBs towards 1 1 stars. 
Equivalent widths are also presented for the following DIBs detected in three sightlines in M 3 1 : /1/14428, 5705, 
5780, 5797, 6203, 6269, 6283, 6379, 6613, 6660, and 6993. All of these M31 DIB carriers reside in clouds at 
radial velocities matching those of interstellar Na i and/or H i. The relationships between DIB equivalent widths 
and reddening (E^_y) are consistent with those observed in the local ISM of the Milky Way. Many of the ob- 
served sightlines show DIB strengths (per unit reddening) which lie at the upper end of the range of Galactic 
values. DIB strengths per unit reddening are found (with 68% confidence), to correlate with the interstellar UV 
radiation field strength. The strongest DIBs are observed where the interstellar UV flux is lowest. The mean 
Spitzer 8/24 //m emission ratio in our three fields is slightly lower than that measured in the Milky Way, but 
we identify no correlation between this ratio and the DIB strengths in M 3 1 . Interstellar oxygen abundances 
derived from the spectra of three M 3 1 H n regions in one of the fields indicate that the average metallicity of 
the ISM in that region is 12 + log[0/H] = 8.54 + 0.18, which is approximately equal to the value in the solar 
neighbourhood. 

Subject headings: dust, extinction - galaxies: individual (M31) - galaxies: ISM - ISM: H n regions - ISM: 
lines and bands - stars: early-type 



1. INTRODUCTION 

Over 400 diffuse inters tellar bands (DIBs) are now known 
dHobbs et al.ll2008t 120091) . but the identity of the carriers has 
remained a mystery since their discovery almost 90 years ago. 
DIBs arise predominantly in the diffuse atomic component of 
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the interstellar medium (ISM) ( ISnow & McCallll2006l) . and it 
is debated whether the y are attribu table to interstellar gas or 
dust (see the review bv lSarr^l2006h . Substructure observed in 
several DIB profiles indicates that the carriers of these DIBs 
are probably large, stable, UV-resistant gas-p hase molecules 
(I Sarre e t al. 1995; E hrenfreund & Fo ing 1996). Likely carrier 
candidates include organic molecules such as polycyclic aro- 
mat ic hydrocarbons (PAHs), fullerenes and carb on chains (see 
e.g. ISalama et al.lll996tlRuiterkamp et al.ll2005l) . 

The study of DIB strengths in extragalactic environments 
allows analysis of the carrier formation and destruction mech- 
anisms under physical and chemical conditions that may be 
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quite distinct from those typically found in the Milky Way 
(MW). Galaxies in the Local Group thus provide us with 
alternative astrochemical laboratories, and opportunities for 
detailed studies of interstellar chemistry over a broad range 
of conditions that may be uncommon or nonexistent in the 
ISM of our Gala xy. Previous ex tragalactic DIB research has 
been reviewed bv lSnowl d2002bl) . Recent research in this area 
has made significant progress using lOm-class telescopes. 
Atoms, molecules and DIBs have been studied in the Small 
and Large Magella nic Clouds (MCs) using optical absorption 
spectroscopy (e.g. lEhrenfreund et alJl2002t |Cox et al.l 120061: 
lCordineril2006l iWeltv et al.ll2006t ICox et al.ll2007h . The ef- 
fects on DIB carrier abundances of the different physical and 
chemical interstellar conditions were examined, including the 
higher gas-to-dust ratios, lower metallicities, different-shaped 
extinction curves and stronger interstellar radiation fields - 
identifying in particular the dependence of DIB strengths (per 
unit reddening) on the interstellar metallicity. Beyond the Lo- 
cal Group, there have been relatively few studies, which have 
genera lly been confined to sightlines probed by bright super- 
novae jRichll 19871; ID' Odorico et al lll989t ISteidel etal.ll 19901 : 
ISollerman et al.ll2005HCox & Patatfeo08l: | Thoene et al.ll2009h 
or background quasars ( e. g. I York et al.l T2006I: lEllison et al.l 
12001 lLawton et al.ll2008l) . 

Desirable attributes of galaxies for detailed optical spectro- 
scopic absorption studies of the ISM include proximity, an 
abundance of OB-type supergiants, Doppler velocities greater 
than ~100 kms~' (to distinguish extragalactic spectroscopic 
features from those arising in the Milky Way foreground), and 
a relatively low foreground extinction. Beyond the Magel- 
lanic Clouds, the obvious candidates for study are the spi- 
ral galaxies M31 and M33. However, these galaxies are 
more than ten times more distant than the MCs, so high- 
resolution spectra are not readily obtained using current tele- 
scopes; surveys of large numbers of stars are only feasi- 
ble at low-to-intermediate spectral resolution. The matter 
in M 3 1 spans a range of radial velocities between around 
and 600 kms~'(see Figure HJ, so the majority of spectral fea- 
tures are easily separable from those in the Milky Way, even 
at modest resolving power. However, as we identify in this 
study, towards the northern edge of M 3 1 the Doppler sepa- 
ration of spectral features can be sufficiently small to cause 
difficulty in separating the MW and M 31 components. 

Recent work using moderate-resolution (Keck DEIMOS) 
multi-object spectroscopy of early-type supe rgiants lead to 
the firs t detections of DIBs in M 3 1 and M 33 (ICordiner et al.l 
l2008allbl) . Unusually large DIB equivalent widths per unit 
reddening (Eb-v), compared to the Milky Way, were found 
in both galaxies. The /I5780 DIB was subsequently found 
to be strong relative to Eb-v in a sample of 17 sight- 
lines in the vicinity o f the M31 OB78 cluster (NGC206) 
(ICox & Cordinerll2008l) . However, these were selected from 
an observed set of 72 sightlines based on their strong DIBs 
relative to the spectral noise, so the results cannot be con- 
sidered to be representative of the typical b ehaviour of DIBs 
in M3 1 . In addition, the blue spectra from ICox & Cordinerl 
(120081) were of very poor quality, preventing the determi- 
nation of precise stellar classifications and reddenings. In 
the present study, emphasis was placed on obtaining higher 
signal-to-noise spectra throughout the optical to facilitate the 
measurement of DIBs and the derivation of stellar spectral 
types. 

In terms of its average gas-to-dust ratio (Nedialk ov et al 
l2000t iBresolin et al]|2002l) and metallicity (e.g. lTrundle et al 



120021) . M 3 1 has more in common with the Milky Way 
than the Magellanic Clouds. However, there is strong evi- 
dence that the ISM in M 31 is different to that of the MW 
and the MCs in the following ways: 1) The M31 star 
formation rate is presentl y about a tenth th at of the MW 
(IWalterbos & Braunl [l994t iKang et all l200l : 2) the inter- 
stellar radiation fiel d is weak in the UV and devoid of far- 
UV radiation (e.g. ICesarskv et al.| [T998t iPagani etailll999l; 
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iMontalto et al.l2009l) - interstellar dust grains have apparently 
not been subjected to the same degree of UV-processing as in 
the MW (Cesa rsky et al.l fl998). which may explain the weak- 
ness of the mid-IR PAH bands in the centre of M 3 1 ; 3) the 
/12175 A UV extinction bump has been measured to be w eak 
and narrow dHutchings et alj|1992t iBianchi et aT] 119961): 4) 
it has a narrow Serkowski polarisation curve (IClavton et al] 
120041) . indicating a different grain-size distribution and/or 
dust-grain composition; and 5) the extinction law is peculiar 
as shown by the anomalously low a verage colour-excess ra- 
tio Eu-b/Eb-v dMassev et al.lll995l) . Accepti ng the curren t 
theories of dust-grain extinction (reviewed by iDrainel 120031) . 
the evidence is consistent with a lack of small graphitic dust 
grains in the ISM of M 3 1 - a conclusion al so reached by 
IXu&Heloul (119941) . iMegier et al.l (120011 120051) identified cor- 
relations between Galactic DIB strengths and the shape of the 
UV extinction curve, which were interpreted as resulting from 
a physical or chemical relationship between the DIB carriers 
and the small grains responsible for extinction. Thus, it is of 
interest to examine whether differences in the size distribu- 
tion of small grains in the M 3 1 ISM compared to the MW 
gi ve rise to any changes in the DIB spectrum. 

iMontalto et al.l (120091) concluded, using UV & mid-IR maps 
of M 3 1 , that the observed dust emission is best represented by 
a model with high PAH abundances (>3% of available carbon) 
and a low UV field. In their scenario, dust is heated predom- 
inantly by an older stellar population, in agreement with the 
low mean intensity of the UV field. The relative lack of UV 
photons in M 3 1 enables examination of the hypothesis that 
a strong UV radiatio n field (on th e level of that found in the 
Milky Way ISM; e.g.\D raine 1978), is required for DIB carrier 
produ ction (see for example. lHerbiglfl995l and iKendall et alJ 
12001 . 

In gener al, DIB equiv alent widths (W) correlate well with 
Eb-v (e-g- lHerbidfl993l) . and may therefore be used to esti- 
mate the amount of dust in a sightline. Thus, DIBs are becom- 
ing important as a measure of reddening in external galaxies 
and D LAs (see for example. lLawton et al. 2008; Ellison et alJ 
120081) . where they may be used to derive the best - and some- 
times the only - estimate of Eb-v- Given the sparsity of 
DIB observation s outside the Milky Way, and the results of 
ICox et alJ (I2006I). IWeltv et al] (I2006I) . ICordiner et al.l J2008al) 
and ICordiner et al.l d2008bl) T"which suggest that extragalactic 
DIBs may deviate from the MW Eb-v relationship, an impor- 
tant goal of our research is to examine the relationship be- 
tween DIB equivalent widths and reddening in M 3 1 . 

In this article we present the largest sample to-date of blue- 
to-red spectra of early-type supergiants in M 3 1 . Following 
the presentation of their spectral classifications and stellar ra- 
dial velocities, we use these data to undertake the first survey 
of DIBs and interstellar Na 1 absorption lines across the disc of 
M 31, and compare the properties and behaviour of the DIBs 
with those measured in the Milky Way. 
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Fig. 1. — Optical image of M 3 1 showing observed fields. Individual tar- 
get locations are plotted in Figure fFTl Credit: Bill Schoening, Vanessa Har- 
vey/REU program/NOAO/AURA/NSF. 



2. TARGET SELECTION AND OBSERVATIONS 

Early-type target stars were selected using UBVRI pho- 
tometry from the Local Group Galaxies Survey (LG GS) cat- 
alogue of 371,782 M31 stars dMassev et al.l 120061) . Photo- 
metric cuts of V < 22, -0.45 < (B - V) < 0.45, and 
—1.2 < Q < -0.4 (where the reddening-free Johnson Q- 
index is Q = (U - B) - 0.72(B - V)), were employed to se- 
lect the brightest OB supergiants from the catalogue. These 
cuts revealed a number of OB associations and star clusters 
across M 3 1 . To check the accuracy of the colour cuts, stars 
in each associatio n were cross-matched with published spec- 
tral classifications (Mass evet al.ll9 86. 1995) and our own un- 
published, low-resolution spectra obtained using the William 
Herschel Telescope. Three fields were selected for observa- 
tion using the Gemini Multi-Object Spectrograph (GMOS) 
on the Gemini North telescope, with the following parame- 
ters optimised: 1) the range of galactocentric distances (i.e. 
relative to the centre of M 3 1); 2) the number of previous OB- 
type classifications in the region; 3) proximity to H n regions 
(the relationship between M 3 1 stellar and H n region metal- 
licities will be examined in a future article). The three fields 
observed are in the fi elds numbered 3,8, and 9 in the survey of 
iMassev et"aH ( 120061) . Our fields 2 and 3 (in Massey's Fields 
8 and 3), include the OB78 and OB8 clusters, respectively 
( I van den Berghll 19641) . The three GMOS fields are shown in 
Figure Q] and cover regions in M 3 1 with galactocentric dis- 
tances of 14.4-18.4 kpc (Field 1), 8.9-9.6 kpc (Field 2), and 
« 5.8 kpc (Field 3). The targets in each field were selected 
to optimise the number of apparently single B-type stars; any 
stars previously classified as binaries, or which looked partic- 
ularly asymmetric in the LGGS images, were excluded from 
the optimisation routine. The orientation of each field was 
optimised to place the majority of stars in the centre to en- 
sure the optimal spectral coverage for each star. Co-ordinates 
and photometry for the final selection of 34 target stars are 
presented in Table Q] 

Observations were carried out between August 2007 and 
January 2008 as part of queue-program GN-2007B-Q-1 16. 



Spectra were obtained in separate blue and red exposures, us- 
ing the B 1200 and R83 1 gratings respectively, to obtain wave- 
length coverage from ~ 4000-7000 A. Using a slit-width of 
1", the mean FWHM of the calibration arc lines was 1.5 A 
across the blue spectral range and 2.3 A in the red. The re- 
sulting resolving power is variable as a function of wavelength 
and ranges from approximately 2500 in the blue to 3500 in the 
red. Seven or eight ~ 30-minute exposures were obtained for 
each field in each wavelength range. Total integration times 
were as follows: Field 1: 3.9 hr (blue), 3.5 hr (red); Field 2: 
4 hr (blue), 4.3 hr (red); Field 3: 4 hr (blue), 3 hr (red). 

The spectra were reduced using the IRAF gemini.gmos 
package! Images were first cleaned of cosmic rays then 
flat-fielded. Individual object images were wavelength- 
calibrated and distortion-corrected using arc lamp exposures, 
then background-subtracted before extraction of the one- 
dimensional spectra. The wavelength calibration accuracy 
was confirmed by measuring the peak wavelengths of the 
Na D sky emission lines. Individual reduced spectra were 
velocity-corrected to the local standard of rest (LSR) frame 
before co-addition. The signal-to-noise per pixel (S/N) of the 
reduced spectra is typically ~ 100 near the centres of the blue 
and red spectral ranges. 

3. RESULTS 
3.1. Stellar classifications 

Previous analyses of luminous B-type supergiants in M 3 1 
found abundan ces comparable to th ose in the solar neigh- 
bourhood (e.g. iTrundle et alj 120021) . The GMOS spectra 
were therefore classified with respect to Galactic standards 
(IWalborn & Fitzpatricklll990t iLennon et all 1 19921) . Spectral 
types and stellar radial velocities are presented in Table Q~l 
those stars with previous spectral classifications are listed for 
comparison with our new classifications in Table [2] Figure [2] 
shows a selection of the blue GMOS spectra, illustrating the 
spectral sequence of the B-type supergiants. 

Reported stellar radial velocities are the mean measure- 
ments of the line-centres of selected H, He and metal absorp- 
tion lines. The velocity for 4029.71 is a notable outlier of the 
Field 2 results and the spectrum shows some evidence for a 
potential second component (e.g. He n /14200, although in the 
opposite velocity sense to what might be expected) - it seems 
likely this is a composite spectrum of more than one star. 

Reddeni ngs (Er y) w ere calculated usi ng the intrinsic 
colours of lWegnerj 1 1994), supplemented by iJohnsonl d 1966b 
for later s pectr al types. The foreground reddening is discussed 
in Sectionl3~4l 

Observational details of the 34 program stars are sum- 
marised in Table Q] We adopt abbreviated target names based 
on the minutes and seconds of right ascension in the LGGS 
catalogue. Thus, stars belonging to each of the three fields are 
distinguished by the first two digits of their abbreviated names 
(39 for Field 1, 40 for Field 2, and 44 for Field 3). 

3.2. DIB sightlines 

Of the stars observed, two are late-type foreground stars 
(4405.89 and 4409.71) and one (3945.82) is a Wolf-Rayet in 
M 31, and therefore unsuitable for our DIB survey. To inves- 
tigate the DIBs in M 31, only the sightlines passing through a 

4 IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the National Science Foun- 
dation. 
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Fig. 2. — Montage of observed M 3 1 supergiants over the range of spectral types 07 Iaf to AO la. Note the trends in the primary diagnostic lines : He i /1447 1 , 
Hen ,U4200, 4541, 4686, Mgn i4481, Sin ii4128, 4130 & Sim AA4552, 4560, 4574. Also labelled in the spectrum of 4029.71 are the Siiv, Nm and Hen 
emission lines. 
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TABLE 1 

Observed target star details 



Field 


Abbr. 


LOGS 


V 


B- V 


Sp. type 


''LSR 


Eb-v 


pM31 

B-V 


1 (F9) 


3936.51 


J003936.5 1 +402235.4 


18.538 


0.049 


B5 la 


-520 ± 2 


0.13 


0.07 ± 0.04 




3944.71 ★ 


J003944.71+402056.2 


18.200 


0.146 


09.7 lb 


-525 ± 5 


0.38 


0.32 ± 0.04 




3945.35 


J003945.35+402115.4 


17.933 


0.325 


B6 la 


-505 ± 5 


0.39 


0.33 ± 0.04 




jy4J.0Z * 


juujy4j.oz-r-4uzjuj.z 




— U.U04 


WIN J + KJ 1 












juujyj4.yj+4Uz4Uj.4 


1 C 77C 
10./ la 


U.U14 


do ia 


^ 1 -i-7 
— J J 1 ± 1 


U.U4 


n no _i_ n c\a 

— U.UZ ± U.U4 




jyjo.y4 


TflfllQ^^ QA _i_/im 1 /IT ^ 
J UU J y J O. y4+ 4UZ 1 4z . 3 


15. /ZO 


o 1 17 
— U. 1 j / 


dU.j la 


—J 10 ± 1 


U.Uo 


U.UU ± U.U4 




■>(lcg 77 


Tnnio^o n i Annua n 


lo.yoo 


u.uyj 


dU. / la 


—J 10 ± 1 




U.z4 ± U.U4 




jyjy.yo 


juujyjy.yo+4Uzz2i . i 


1 Q £7^ 

lo.OZj 


n in 
— u. 11/ 


7 TK 

\Jy. 1 ID 


—J 14 ± 4 


n 1 1 
U. 1 1 


U.Uj ± U.U4 


Z ( vJjd la) 


4UZO. / 4 


Tn0409^ 74-4-40441 8 8 
JUU'tUZO. / 4+4U44 lo.o 


1 8 387 
1 O.JO / 


— u.u / 1 


DJ Id 


S7Q -t- 7 

J /y ± 1 




(in'; 4.(1 HA 

—U.UJ ± U.Ut 




4027 18 


T004097 1 R-t-ATldW, ? 


18 009 


—0 050 


B3 Ia 


—571 + 2 


03 


— rn + 04 

_r_ \J.\J^ 




4029 57 


TftfWnQ S7+4044S7 s 


18 192 


w.uzu 


B8 Ia 


_SQO -t- 7 


01 


o ns -1- o 04 

— U.UJ _£ U.Ut 




4029 71 


T00409Q 7I-+-40449Q 8 


18 561 


—0 229 


07-7 5 Iaf 


—651 + 5 


08 


09 + 04 




403CU8 * 


J004030.28+404233.1 


17.357 


-0.042 


Bl Ia 


-579 ± 3 


0.15 


0.09 ± 0.04 




4030.32 * 


J004030.32+404404.8 


18.950 


-0.063 


B5 Ia 


-578 ± 3 


0.02 


-0.04 ± 0.04 




4030.84 


J004030.84+404348.6 


18.827 


-0.093 


B2 lab 


-566 ± 3 


0.07 


0.01 ± 0.05 




4030.94 


J004030.94+404246.9 


18.866 


-0.147 


09.5 lb 


-546 ± 5 


0.09 


0.09 ± 0.05 




4031.52 * 


J004031.52+404501.9 


18.923 


-0.147 


B0.5 Ia 


-591 ±4 


0.05 


0.03 ± 0.05 




4032.28' ★ 


J004032.28+404506.7 


17.916 


-0.001 


B8Ia 


-578 ± 4 


0.03 


-0.03 ± 0.04 




4032.88 * 


J004032.88+404509.9 


17.952 


-0.103 


B2.5e 


-590 ± 8 


0.04 


-0.02 ± 0.04 




4032.92 


J004032.92+404257.7 


18.193 


-0.092 


B0.5 Ia 


-594 ± 2 


0.11 


0.05 ± 0.05 




4034.61 


J004034.61+404326.1 


18.669 


0.150 


Bl Ia 


-561 ± 2 


0.34 


0.28 ± 0.04 




4034.68 


J004034.68+404509.9 


18.959 


0.037 


AOIa 


-582 ± 3 


0.03 


-0.03 ± 0.04 




4034.80 


J004034.80+404533.5 


18.378 


0.254 


B8Ia 


-584 ± 2 


0.28 


0.22 ± 0.05 




4037.92 


J004037.92+404333.3 


18.657 


0.056 


B1.5Ia 


-571 ±3 


0.23 


0.17 ±0.04 


3 (OB8) 


4405.89 f 


J004405.89+413016.3 


17.309 


0.847 


K0 




0.03 






4406.87' 


J004406.87+413152.1 


19.019 


-0.116 


B1.5Ia 


-59 ±4 


0.05 


-0.01 ± 0.04 




4408.04 


J004408.04+413258.7 


18.708 


-0.017 


B3Ia 


-105 ± 3 


0.11 


0.05 ± 0.04 




4408.36 


J004408.36+413210.2 


19.339 


-0.098 


B0.7 Ia 


-104 ±4 


0.10 


0.04 ± 0.05 




4408.79 * 


J004408.79+413142.0 


19.602 


-0.050 


08-9 II 


-145 ± 1 


0.24 


0.18 ±0.06 




4408.94 


J004408.94+413201.2 


18.356 


-0.088 


B1.5Ia 


-127 ±5 


0.08 


0.02 ± 0.04 




4409.52 


J004409.52+413358.9 


19.302 


-0.068 


B2Ib 


-97 ±4 


0.12 


0.06 ± 0.05 




4409.7 l f 


J004409.71+413247.2 


16.834 


0.547 


F8 




0.00 






4412.17 


J004412.17+413324.2 


17.330 


0.345 


B2.5 Ia 


-110 ± 5 


0.49 


0.43 ± 0.04 




4417.80 * 


J004417.80+413408.0 


19.397 


-0.055 


09.5 lb 


-94 ± 1 1 


0.19 


0.13 ±0.05 



Note. — Local Group Galaxies Survey (LGGS) designations and photometry are from|Massev et al. (2006), abbreviated 
names used in this paper (Abbr.) are also given. Total line-of-sight reddenings (Eb-v in magnitudes) and radial velocities (vlsr 

kms" 1 ) are given. M31 reddenings (E™^) assume a foreground reddening of 0.06 mag as discussed in Section l3~4l Stars 
suffering from possible spe ctral and/or photom etric contamination due to blending with nearby objects (as determined from 
the LGGS V-band images; Massey et al. 2006), are labelled with a star. The spectrum of 3945.82 was kindly classified by 
Prof. P. Crowther (private communication). Two stars in Field 2 have negative E^_y values outside of the error bars, which may 
be indicative of a problem with their photometry, perhaps due to unseen multiplicity or variations in the assumed foreground 
reddening. 

* Possibly blended with other objects. 
' Stars used as telluric standards. 
f Foreground stars. 
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TABLE 2 

New classifications compared to previous work 



Star 


Alias 


Sp. Type 


Previous 


3945.82 


OB 135-1 


WN5 + 07 


WN 2 


4029.71 


OB78-231 


07-7.5 Iaf 


08.5 1(f) 1 


4030.28 


OB78-277 


Bl la 


B 1.5 la 3 


4032.88 


OB78-478 


B2.5e 


B1.5 la 4 


4032.92 


OB78-485 


B0.5 la 


B0 I 5 


4037.92 




B1.5 la 


B0 1 6 


4406.87 


OB8-7 


B1.5 la 


Bl I 1 


4408.36 


OB8-25 


B0.7 la 


Bl I 1 


4408.94 


OB8-34 


B1.5 la 


Bl la 1 


4409.52 


OB 10-43 


B2 lb 


B8Ia' 


4412.17 




B2.5 la 


B5: I + Hn 6 


4417.80 


OB9-176 


09.5 lb 


Early B 1 


References. — (11 IMassev et al.1 


11995): (2) 



I Massev et al] (l9EB i: ( 3) IBianchi et al.1 



1990); (6) 



Trundle et al. 12002); (5) Humphreys et al 
Massey et al. i200g). 



significant quantity of ISM were selected for further analysis. 
Each spectrum was exa mined for the presen ce of DIBs, by 
reference to the survey of lHobbs et al.l (|2008). No DIBs were 
detected towards stars with a foreground-corrected reddening 
of < 0.05 mag. This is as expected given the spectral 
noise, and the assumption that M 3 1 DIB strengths are similar 
to those typically observed in the Milky Way: For the typi- 
cal signal-to-noise (S/N) of the spectra near 5780 and 6283 A 
(100 and 80, respectively), using the ap proximate minimum 
DIB W detection threshold equation of iHobbs et al] d2008l) 
Wma = l-064xFWHM/(S/N)), yields respective lower de- 
tection limits of 22 and 63 mA for these DIBs. Scaled to a red- 
dening of Eb-v = 0.04, the mean Galactic /L/15780 and 6283 
DIB W/E R y values, ar e 21 and 47 mA respectively (using 
data fro m iHerbid [19931 iThorburn etail 120031 iMegier et al] 
120051 and ICordinerll2006l) . which is consistent with their non- 
detection in our spectra at M31 reddenings of 0.04 and less. 
The sample of 14 stars with E^l > 0.05 is listed in Table 
H with the addition of three targets (3956.94, 4032.28 and 
4406.87), which were used as telluric standards for Fields 1, 
2 and 3, respectively. 

3.3. Sodium D lines: column densities & radial velocities 

At the resolution of the GMOS spectra (100 kms -1 in this 
wavelength region), the interstellar Na D absorption lines 
in M31 and the MW can be modelled using single Gaus- 
sian interstellar cloud components. Both lines of the Na D 
doublet were fitted simultaneously using the vapid routine 
dHowarth et alj|2002l) . to produce models of the Galactic and 
M31 Nai distributions in each sightline. Fitting both lines 
of the sodium doublet simultaneously helps to constrain the 
Doppler b parameter in data at this relatively low spectral res- 
olution. The resulting interstellar Na i cloud radial velocities, 
Doppler b parameters and column densities are summarised 
in Table [3] Monte Carlo error estimates were generated by 
refitting each cloud model to 100 different replicated spectra, 
each with random Poisson noise added, and taking the +68th 
percentiles of the resulting (free) parameter ranges. 

The Doppler shift of Field 3 is comparable to the instru- 
mental resolving power, and as a result, the MW and M 3 1 
absorption components are blended which hinders the accu- 
rate derivation of the individual component parameters. To 
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Fig. 3. — Observed continuum-normalised Nai spectra (histograms) with 
best-fitting interstellar cloud models overlaid. The velocity scale at the top 
is relative to the Na Di rest wavelength. Stellar He ,15876 lines have been 
removed from the observed spectra using Gaussian fits. Milky Way (MW) 
and M 3 1 absorption components are marked. Dotted lines show the Na D 
rest wavelengths. 'Denotes stars used as telluric standards. 
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TABLE 3 

Best-fit interstellar Na i cloud model parameters for the Galactic (MW) and M 3 1 ISM 







MW 






M 31 




Sightline 


V 


b 


N 


V 


b 


N 


3936.51 
3944.71 
3945.35 
3956.94' 
3958.22 
3959.98 


-0 4 +41 

"■ -5.6 

11 5+ 5 - 8 

-6.2 

-4.55+61 
-46+« 

-45.6+.™ 

-41 2+ 4 - 8 

*'■ -6.5 


4 8+ 37 
10.4+ 37 / 

4Q 4+11.9 
^"-13.9 

3 - 8 -a9 

5 2+ 24 

J - -2.1 

5 9+ 3 9 


12 43+ 028 
-0.22 

12 14+° 18 

12 09+° 03 

12 77+ - 45 
-0.27 

19 7Q+0.69 
lz -' y -0.26 

19 47+ - 27 
1Z - -0.18 


-526.0+| 4 

-518.5+1" 

-514.0+|-| 

-521.4+ 7 ; 8 

-516 6+ 4 - 1 
jiu.u _ 46 

-502.5+. 3 ; 7 


7 g+i.2 
'•'-1.1 

10 5+ L6 

11 l+'o 

"■'-0.9 

3 7+ 3 - 1 

7 0+ 22 
'•"-1.4 

3.0+.0- 4 


12.90*» 
13 57+0- 38 

■ -0.29 

13 16+0" 

1J '"-0.10 

12 36+0- 56 

1 3 9Q+0-45 
1J - Z '-0.41 

13 56+0" 63 
U,J -0.42 


4030.28 
4032.28' 
4034.61 
4034.80 
4037.92 


-2 0+ 30 

z - -3.8 

-34 l+ 3 -° 

J -3.2 

-35 7+ 5 " 4 
JJ -'-6.3 

-49.6+3;5 
-29.8+ 4 ;^ 


51 3+ 7 - 3 
ji.j_ 70 

10 2+ 53 
'"-2.8 

o c+21.1 

8 - S -6.3 

26 0+ 8 - 9 
z "'"-10.2 

ao 7+13.8 

^-15.2 


1 9 aq+0.03 
1Z ' J -0.02 

12.58+°;» 
12 41+° 24 

1Z -<+'-0.19 

12 50+ 04 
12 31+0-04 

lz - J '-0.04 


-578.5+j 8 

-57 1 7+ 22 - 5 
' '-'-16.3 

-552.2+| 8 
-548.S+ 2 - 8 
-554 2+ 3 - 1 

JJH - -3.6 


46.8+" 
58.2+^ 

62-21:? 
11.6+g; 6 
7 4+2.0 

'■ -1.6 


1 9 94+0-03 

1 1 79+ 012 
"■"-o.io 

1 9 7S+0.02 
lz - /S -0.02 

13.08+0- 22 
13.26:0|i 


4406.87' 
4408.04 
4408.79 
4409.52 
4412.17 
4417.80 


-34.0+ 3 ! 
-32 0+ 5 - 2 

JZ -"-4.2 
4 4+10.2 
H " -11.6 

-23.2 f 

-23.2 f 

-40A +6 ,°. 
-6.4 


16 9+ 75 

'"■'-6.0 
99 7+10.2 
■"■•'-11.0 

10 5+ 15 5 

'"-5.6 

19.9 f 
19.9 f 

93 4+9.4 
ZJ - -6.0 


19 34+0.04 

12 55+ 012 

-0.05 

12.81+^ 
12 43+0- 05 
12.28+0- 2 o 
12 60+0 .12 


-121.0* 
-151.3+ 2 -0 
-90.5+10- 7 
-119.2+11 

-122.9+0- 7 


12.0 f 
10.9+.*?, 
7 7+ 3 - 2 

'■'-3.5 

16 2+ 2 - 6 

'"■ -3.0 
18.5+0-9 

7 5+ 16 
'■ J -2.4 


< 11.4 

13.48+Of 
13 10+0- 95 

■ -0.36 

13 16+ 018 
13 25+0- 05 
13 22+0- 47 



Note. — Radial LSR velocities (v) and Doppler b parameters are in km s . Column 
densities (N) are in log cm - -. Statistical Monte Carlo errors are given. The MW model 
parameters may be uncertain due to contamination of the Na D lines near v = by sky 
subtraction residuals. M31 column densities should generally be considered as lower limits 
due to the likely presence of unresolved saturated structure in the Na D lines. 

f Values held fixed during the fitting. 
' Telluric standard stars. 



obtain physically realistic results for 4409.52 and 4412.17, 
the velocity and Doppler b parameter of the MW component 
had to be held fixed during the fitting (at v = -23.3 kms~' 
and b = 19.9 kms~', which are the averages for the MW 
components towards the other stars in Table [3j. Towards the 
telluric standard 4406.87, there was no detectable M 31 com- 
ponent of Na i; the only absorption detected is at a mean radial 
velocity of -32.0^ 2 kms~' which is unlikely to be in M 31 
and is consistent (in velocity and column density) with the 
other MW foreground components observed towards the stars 
in Fields 1 and 2. An upper limit on the M31 Nai column 
density of 2.8 x 10 11 cirT 2 was derived for this sightline as- 
suming the radial velocity and Doppler width of the gas to be 
the average of those values measured towards the other stars in 
Field 3. It must be noted, however, that all of the Field 3 Na i 
fit parameters are subject to additional uncertainty due to the 
likely contamination of the spectra by sky-subtraction resid- 
uals. Strong telluric Nai emission features dominate these 
GMOS spectra at velocities +100 kms -1 with respect to the 
D-line centres, and their imperfect subtraction during data re- 
duction thus leads to potential uncertainties in the spectra over 
this velocity range. 

The observed sightlines pass through the halo of M31 
and are expected to intercept multiple interstellar clouds with 
component structure that is unresolved in our spectra. Thus, 
the derived (single-component) interstellar cloud models pro- 
vide no information about the detailed kinematics of the M 3 1 
ISM except for the overall velocity width and peak radial ve- 



locity of the Nai gas distributions. The column densities de- 
rived from these broad, single-component fits should gener- 
ally be considered as lower limits due to the likely presence 
of unresolved saturated structure in the Na D lines. 

Since sodium lines can be present in the background stellar 
spectra, we investigated the possible stellar Na D contribu- 
tions to our observed spectra. Theoretical spectra for a range 
of stellar types spanning those o bser ved were calcu lated using 
SYNSPEC dHubenv & La nz 2000) and|KurucJ dl993|) model at- 
mospheres with solar metallicity. Stellar contributions to the 
Na D line equivalent widths were found to be negligible com- 
pared with the interstellar component for all the stars listed 
in Table [3j only in cooler stars (of type A and later) does the 
stellar contribution become important. 

The measured M31 interstellar sodium radial velocities 
match the velocities of nearby M 3 1 stars and planetary neb- 
ulae, as can be seen in Fig. |U This demonstrates with little 
ambiguity that these clouds are located in the main disc of 
M31. 

3.4. Foreground reddening and DIBs in the Galactic halo 

The mean foreground Na i column d ensity towards our sur- 
veyed M31 stars is 3.2 x 10 12 cm -2 . ISembach et"aT1 (1 1993b 
measured the Na i column density in lines of sight through the 
Galactic halo, probing gas that should be similar to the MW 
foreground ISM in our sightlines. The correlation between 
Nai and Eb-v permits an estimate of the f oreground red- 
dening. Performing a least-squares fit to the ISembach et alj 
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Fig. 4. — Kinematic (position- velocity) map of M 3 1 along the galaxy's ma- 
jor axis, sho wing the emission-lin e objects (mostly planetary nebulae) from 
the survey by 'Merrett et al. (2006) and our measured interstellar Nai veloci- 
ties. Data points for each GMOS field are plotted using a different symbol. 
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( 1 19931) data yields a linear relationship: E B -v = (N(Na i) + 
1.41 x 10 n )/3.79x 10 13 , from which we derive a foreground 
reddening of Eb-v = 0.088 mag. However, there is consider- 
able scatter in the Nai vs. Eb-v data - the mean RMS error 
on Eb-v is 0.053 - and there is likely to be contamination of 
our Na i profiles near v = due to sky-line subtraction resid- 
uals. These factors makes the derivation of foreground red- 
dening from Na i qui t e uncertain. From the LAB H i survej0 
dKalberla et al.ll2005t iHartmann & Burtori [19971) . the spectra 
obtained at the points nearest the centres of our three fields 
give foreground Hi column densities of 5.75 x 
(Field 1), 6.07 x 10 20 cnr 2 (Field 2), and 6.68 x 
(Field 3), assuming optically thin emission. Using Equa- 
tion 7 of iBurstein & Heilesl (1 1978b (which assumes a con- 
stant gas-to-dust ratio) to calculate the foreground redden- 
ing, the average value is Eb-v = 0.068 mag. The modi- 
fied lBurstein & Heilesl relation utilised by the NED database^ 
results in Eb-v =0.0625 mag, which is in good agreement 
with the foreground reddening towards M31 estimated by 
ISchlegel et al.l (fl998) from the median Galactic dust emission 
in the surrounding regions (Eb-v = 0.062 mag). We there- 
fore adopt a MW foreground reddening of 0.06 + 0.02, which 
is used in the calculation of the foreground-corrected redden- 
ings (Eg 3 y) given in Table Q] Quoted errors on E^Jy include 
an uncertainty of 0.02 in the foreground reddening, plus the 
uncertainty on the intrinsic stellar B - V values (equal to the 
uncertainty of a single spectral subclass). 

As shown in Figure [5] the i5780 DIB is detected in 
several of the observed sightlines at a wavelength match- 
ing the Galactic /t5780 rest wavelength of 5780.37 A (see 
iGalazutdinov et ai1l2000l for example). This absorption fea- 
ture most likely originates in the Galactic halo. The average 
foreground /15780 equivalent width is 40 + 10 mA, as mea- 
sured from the co-added spectrum of all the stars from Field 
1 and 2 listed in Table [1] (excluding the Wolf-Rayet and those 
sightlines in Field 3 with blended MW and M 3 1 Nai cloud 
components). Using a Gaussian fit, the mean radial velocity 
of the i5780 carrier-gas was found to be 0.5 + 7.3 kms . 

5 http://www.astro.uni-bonn.de/~webrai/english/tools_labsurvey.php 

6 The NASA/IPAC Extragalactic Database (NED) is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, under contract 
with the National Aeronautics and Space Administration. 



The mean Galactic /15780 equivalent width per unit redden- 
ing is [W(5780)/EB-v]mean = 527 mAmag -1 , with a stan- 
dard deviation of 230 mAmag -1 (see Section [3~2l for details 
of the literature data used; the outliers towards p Oph and 
the Orion nebula have been excluded). Using this ratio, the 
average Milky Way foreground /15780 equivalent width mea- 
sured towards M 3 1 corresponds to a foreground reddening of 
Eb-v = 0.08^^4 mag, which is consistent with the value of 
0.06 adopted in this study. 

3.5. M 31 diffuse interstellar bands: equivalent widths and 
radial velocities 

We searched the spectra of each of the sightlines with 
foreground-correcte d reddening E¥ 3 v j > 0.05 for the DIBs 
listed in the survey o flHobbs et ail (120081) . For the ,Lt6283 and 
6993 DIBs, which are partially obscured by overlapping tel- 
luric lines, corrections were performed by division with spec- 
tra from lines-of-sight with negligible Eg_y selected from the 
same field (designated with superscript t in Table [TJ. The 
AA4428, 5705, 5780, 5797, 6203, 6269, 6283, 6379, 6613, 
6660 and 6993 DIBs were detected in at least one sight-line 
each, with Doppler shifts matching the M 31 N a i and H i ve- 
locities (Hi data are presented in Section [3~6b . The /15849 
and /16445 DIBs were tentatively detected in the spectrum of 
3944.71 only. The AA4501, 4726, 6196, 6376 and 6532 DIBs 
were too weak to be detected in any of the spectra, which 
is consistent with the assumption of MW-like strengths (per 
unit reddening), given the S/N and low resolving power of 
our spectra. Severe contamination by telluric absorption pre- 
vented analysis of the AA6886, 6919, 7224 and 7334 DIBs. 
Equivalent widths were measured by integration across a lin- 
ear or quadratic fitted continuum. Errors on DIB equivalent 
widths measured this way were determined by crAA, where 
<x is the RMS noise in the continuum and AA is the full- 
width half-maximum of each DIB. As a result of the rela- 
tively low signal-to-noise and low dispersion of the spectra 
compared with othe r contemporary Galactic DIB studies (e.g. 
iMcCall et al.ll2.010b - additional uncertainties in the measured 
DIB equivalent widths arise due to the difficulties in establish- 
ing the true full-widths of the DIBs over which to perform the 
integration. For example, in Figure the presence of over- 
lapping weak stellar/noise features at the wings of the /L5780 
DIB observed towards 3944.71, 4034.80 and 4412.17, cause 
depressions in the continua that cannot be distinguished from 
the DIB by eye, and thus artificially boost the measured equiv- 
alent widths. Conversely, the broad, shallow, Lorentzian-like 
wings of the /16283 DIB (shown in Figure [6]), are lost in the 
noise, resulting in likely underestimation of the full-width of 
this feature. Therefore, to provide a second, more accurate 
measurement of the AA5180 and 6283 radial velocities and 
equivalent widths that takes the complex profiles of the DIBs 
into account, Galactic DIB template spectra were fitted to the 
observed M 3 1 spectra. 

The Galactic DIB templates were derived from high resolu- 
tion (R = 58, 00), high S/N (~ 2000) spectra of B l Sco (see 
Cordiner 2006), shifted to the interstellar rest frame. These 
were convolved with a Gaussian interstellar cloud model, then 
with the GMOS spectral PSF (assumed Gaussian), and re- 
binned to the dispersion of the GMOS spectra. The interstellar 
cloud model was optimised using the same n onlinear least- 
squar es algorithm employed by the vapid code dHowarth et al] 
120021) . This technique has the benefit that the statistical un- 
certainties in the measurements of the equivalent widths and 
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TABLE 4 

,15780 and ,16283 DIB equivalent widths and velocities, M31 Hi21 cm column densities, gas-to-dust ratios, interstellar NUV radiation field strengths and Spitzer 8/24 /jm flux 

RATIOS FOR OBSERVED SIGHTLINES 



W(5780) W(6283) 



Sightline 


Fit 
mA 


Int. 
111 A 


Fit 
m A 


Int. 
m A 


v(5780) 

, -1 
kms 


v(6283) 

1 -1 

km s 


N(H 1) 
10 cm 


G 


Inuv 
Draine 


8/24 11m 
flux 


W(5780)/£j,_ y 

A -1 

A mag 


W(6283)/£ J 

* -1 
A mag 


3936.51 
3944.71 
3945.35 
3958.22 
3959.98 


04+33 

°^-32 

318^5 
"202+ 22 
21 1+ 33 

^"-48 

< 29 


89 (19) 
359 (17) 
b 302 (14) 
207 (35) 
<42 


979+96 

m+J* 

587!™ 
439+98 

1 1 ' -60 


217(55) 
544 (48) 
414(52) 
419 (60) 
124 (54) 


-5 1 7+ 36 s 

Jlu " -36.0 

-512 4 +7 5 
-515 4+ 12 - 2 

J1J '^-11.7 

-517 9 +ls l 

-16.7 


- 498 - 5 -49J 

-524.8+I 3 - 6 
-510 6 +15 6 
-511.9+ 2 " 
-481.3+|; : | 


2.8 
4.5 
4.2 
4.0 
2.6 


20 
7.1 
6.4 
8.5 
25 


0.27 
0.48 
0.35 
0.35 
0.59 


1.16 
0.63 
0.94 
1.78 
0.71 


6 +0 - 2 
9+ 0J 

u - -0.3 

< 2.2 


4.0+H 

■j 1+O.6 
Z ' A -0.4 

1 9+ - 9 
-0.6 

<) 9+13.8 
i - x -1.6 


4030.28 
4034.61 
4034.80 
4037.92 


<41 

215^ 
< 39 


< 28 
165 (25) 
238 (22) 

<49 


<42 
3 75 +i05 

383 + ™ 


< 112 
405 (65) 
326 (54) 
123 (68) 


-514.9! 333 
-584.511" 


-601 3 +310 

oul -30.4 

-543 5 +25J 
-600.1^;, 


2.8 
2.3 
1.3 
3.3 


16 
4.0 
2.9 

10 


3.1 
2.9 
1.9 
0.97 


0.90 
1.20 
1.15 
1.38 


< 0.9 

4 +0 - 2 

1 0* 0A 
<0.3 


< 0.9 
1 3 +0 - 7 

'■ J -0.4 

1 7 +0 - 9 
'■ -0.6 

1 l +u 

'■'-Of) 


4408.04 
4408.79 
4409.52 
4412.17 
4417.80 


< 30 

b 145+74 

< 38 
190^ 

79 -41 


<40 

< 88 

< 57 
256 (10) 

< 57 


< 243 
189+]™ 

<65 
491+83 

194+>f 


< 142 
171 (104) 

< 168 
423 (53) 
221 (76) 


-88.lt 2 

_99 3+57.8 
,7 -''-48.6 


-93.8!f 7 ;S 

-82.8^ 
-107.0!«; 3 


1.6 
1.2 
1.4 
2.0 
1.7 


15 
3.4 
1 1 

2.3 
6.6 


0.79 
0.83 
0.59 
1.1 
0.48 


1.16 

1.05 
1.26 
0.90 
0.76 


< 2.3 

< 3.2 


< 2.7 

1 r" 
'■'-0.7 

< 5.4 

1 i +L9 
'■'-0.7 

1 6 +3 -° 
'■"-l.o 



Note. — Fitted and conlinuum-inlegraled DIB equivalent widlhs arc given, labelled Fit and Int.. respectively. The errors {+) on the integrated equivalent widlhs are crAA and are given in brackets, with upper limits calculated as 2itA,( . Errors on the 
gas-to-dust ratios (G - \ N(Hi)/E^^ ) are ±100% due to probable contamination by background Hi (see Sectior <3.6| . 

a Units of G are 10 21 cm -2 mag -1 

k DIB equivalent widths probably affected by line blending. 



velocities of weak DIBs are minimised, but has the drawback 
that the intrinsic M 3 1 DIB profiles and rest wavelengths are 
assumed to closely match those of the Milky Way DIB stan- 
dard star /3 l Sco. If this were not the case, unknown system- 
atic errors could occur, but there is no evidence to suggest 
that M 3 1 DIB profiles differ from those of /i 1 Sco. Profiles of 
some Galacti c DIBs have been shown t o exhibit variations in 
substructure dGalazutdinov et ail 120021) . So far, however, all 
observed variations have been small and are negligible com- 
pared to the size of the GMOS resolution element and the 
level of noise in our spe ctra. The reader is referred to Section 
3.5.2 of lCordineri ( l2006l) for a more detailed discussion of the 
merits of the DIB fitting technique. 

The fitted and integrated equivalent widths are shown in Ta- 
ble[5] and match each other to within the error bars. However, 
for four out of seven cases (for which the DIBs are not affected 
by obvious blends, denoted by superscript b), the fitted values 
of W(5780) are systematically less than the integrated values. 
For W(6283), the fitted values are generally greater than the 
integrated values. These differences are as expected given the 
increased accuracy afforded by using known DIB profiles to 
better constrain their full-widths, as explained above. 

Error estimates for the /15780 and /I6283 equivalent widths 
and velocities were derived by the same Monte Carlo tech- 
nique as used for the Nai spectral modelling; 500 different 
replicated spectra were produced for each DIB with (random) 
Poisson noise added to each, the cr-value for which was taken 
to be equal to the RMS noise of the local continuum. Errors 
due to the uncertainties inherent in the spectral continuum rec- 
tification were accounted for by multiplicatively scaling each 
replicated spectrum by a random factor of between 1 — cr and 
1 + cr. The DIB model parameters were refitted for each repli- 
cated spectrum and individual parameter errors were obtained 



from the +68th percentiles of the respective parameter ranges. 

As can be seen from the data in Table [5] the error estimates 
on the DIB equivalent widths derived from Monte Carlo fit- 
ting (in the column labeled 'Fit.') are typically larger (by 
~ 50%) than those derived using crAA (column labeled 'Int.'). 
Although our Monte Carlo errors do not take account of in- 
terloping stellar features (these are considered separately; see 
below), they are expected to be more reliable than the crAA 
errors, which assume a somewhat arbit rary + 1 <x erro r on th e 
continuum placement, as described by iMcCall et al.l d2010l) . 
and do not take into account the non-Gaussianity of the DIB 
profiles, which is especially important for /16283 with its very 
broad, Lorentzian-like wings. 

The spectral regions spanning the /15780 and A5791 DIBs 
are shown in Figure [5] for the fourteen stars in our M 3 1 red- 
dened sample (plus three telluric standards). The spectral re- 
gions covering /16269 and /16283 are shown in Figure|6](these 
spectra have been divided by the spectra of the telluric stan- 
dards from each respective field). M 31 components of /15780 
are detected in 10 sightlines and /t6283 in 11 sightlines. The 
MW foreground /16283 components are effectively removed 
during telluric division because of their presence in the tel- 
luric standard spectra. 

A total of 1 1 different DIBs (^4428, 5705, 5780, 5797, 
6203, 6269, 6283, 6379, 6613, 6660, and 6993) were detected 
in M31 towards the three most reddened targets 3944.71, 
3945.35 and 4412.17 (data presented in Tableland Figs.g] 
to [TBI. Figu res [8] to [13] also show reference MW DIB spectra 
(from ICordineri 120061) . as well as telluric standard spectra to 
illustrate the likely level of contamination by stellar and tel- 
luric features. 

To account for possible additional uncertainties in the DIB 
equivalent widths given in Table [5] due to the presence of 
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Fig. 5. — Observed continuum-normalised spectra of the /15780 region (his- 
tograms). /!5780 components arising in MW and M 3 1 have been fitted using 
the Galactic sightline /?' Sco as a template, (fitted profiles overlaid in red). 
Velocity scale at the top is relative to the Galactic /15780 rest wavelength 
of 5780.37 A IGalazutdinov et al.ll200(]) (plotted with a vertical dotted line). 
The mean M 3 1 interstellar Na i absorption wavelengths are plotted with ver- 
tical dashed lines. For the Field 3 stars, the peak H i velocities are also shown 
(dot-dashed lines). 'Denotes telluric standard stars. 
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Fig. 6. — Observed continuum-normalised, telluric-corrected spectra of the 
/i6283 region (histograms), with fitted DIB profiles overlaid (red curves). Ve- 
locity scale at the top is relative to the Galactic /16283 rest wavelength of 
6283.85 A IGalazutdinov et al. 2000) (plotted with a vertical dotted line). 
The mean M 3 1 interstellar Na i absorption wavelengths are plotted with ver- 
tical dashed lines. 
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Fig. 7. — Observed continuum-normalised spectra of the A44 2S region (his - 
tograms). Upper panel shows the intrinsic MW /14428 profile I Snow 2002a). 
Velocity scale at the top is relative to the Galactic DIB rest wavelength of 
4428.39 A (dotted line). M31 Nai velocities are plotted with dashed lines. 
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Fig. 9. — Observed continuum-normalised spectra of the /16203 region. Up- 
per p anel shows MW comparison DIB spectrum towards /?' Sco jCordinerl 
2006J). The DIB rest velocity is shown with a dotted line and the M31 Nai 
velocities are plotted with dashed lines. For 4412. 17, the mean H i velocity is 
shown with a dot-dashed line. 
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Fig. 8. — Observed continuum-normalised spectra of the /15705 re gion. Up- 
per p anel shows MW comparison DIB spectrum towards p Oph A (Cordineil 
I2006D . The second panel from the top shows the spectrum of one of the 
tellu ric sta ndard stars. Velocity scales at the top of this Figure and Fig- 
ures |9|13| ar e relative to the Galactic DIB rest wavelengths published by 
Galazutdinov et al. (2000). The mean M 31 interstellar Nai absorption wave- 
lengths are plotted with vertical dashed lines. The spectrum of 4412.17 is 
dominated by stellar features at this wavelength. 
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Fig. 10. — Observed continuum-normalised spectra of the /16379 region. 
Uppe r panel shows MW comparison DIB spectrum towards /?' Sco ICordinerl 
2Q06). The DIB rest velocity is shown (dotted line), and the M31 Nai ve- 
locities are plotted with dashed lines. For 4412.17, the mean Hi velocity is 
shown with a dot-dashed line. The strong absorption features around 6360 
and 6390 A are due to stellar S n and Ne i. 
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Fig. 11. — Observed continuum-normalised spectra of the ,16613 region. 
Uppe r panel shows MW comparison DIB spectrum towards p 1 Sco ICordinej 
12009) . The DIB rest velocity is plotted (dotted line), and the M31 Nai ve- 
locities are plotted with dashed lines. For 4412.17, the mean Hi velocity is 
shown with a dot-dashed line. 
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Fig. 13. — Observed continuum-normalised spectra of the /16993 re- 
gion. Upper pa nel shows MW comparison DIB spectrum towards p Oph A 
jCordinej|2006ft . The DIB rest velocity and M 31 Nai velocities are plotted 
with dotted lines. 
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Fig. 12. — Observed continuum-normalised spectra of the /J6660 region (his- 
to grams). Upper p anel shows MW comparison DIB spectrum towards p Oph 
A )Cordineil2006l) . The DIB rest velocity and M 3 1 Na i velocities are plotted 
with dotted lines. The strong stellar line at the red end is due to He i. 



overlapping stellar lines, we have estimated the possible stel- 
lar line contributions. The table shows the total equivalent 
widths of stellar features that might be overlapping the DIBs, 
calculated from synthetic B0 la and B9 la stellar spectra (at 
temperatures of 28,000 and 1 0,000 K, respectively), using the 
method described in Section 13.31 The line list used in the 
spectral modeling includes transitions of the following ele- 
ments and their ions: H, He, C, N, O, Na, Mg, Si, Ca, Cr, Mn 
and Fe. The calculated stellar line contributions to measured 
DIB equivalent widths are generally small (except for the very 
broad /14428 DIB which is contaminated by the O n blend near 
4415 A), and in most cases are negligible compared with the 
errors arising from continuum placement and statistical noise. 

As a result of the close blending between the MW and M 3 1 
interstellar gas components due to t he lo w Dopp ler shift of 
Field 3 (discussed further in Sections [3.3| and [3~6l l, there may 
be some ambiguity as to in which galaxy the Field 3 DIBs 
originate. The main discriminator we use to identify the loca- 
tion of the DIB carriers is their radial velocity compared with 
the Na i component structure. The separation of the MW and 
M31 DIB velocities is clear in all cases apart from 4412.17 
and 4417.80, for which the DIB velocities appear to be inter- 
mediate between the M31 and MW Nai components. How- 
ever, a number of factors lead us to conclude that the DIBs 
observed in these sightlines originate in M 3 1 : 1 ) The DI B ve - 
locities match the peak M31 Hi velocities (see Section [3~6b . 
whereas we detect a negligible contribution to the DIB equiv- 
alent widths at radial velocities corresponding to H i in the 
MW (at v/^ < 45 kms 1 ; see iBraun et"aT1 120091): 2) LAB 
survey data towards the targets in question dKalberla et alj 
l2005UHartmann & Burtonl 19971) show a relatively small fore- 
ground Hi component, such that the W(5780)/W(Hi) ratios 
towards 4412.17 and 4417.80 would be very large if this 
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TABLE 5 

DIB EQUIVALENT WIDTHS FOR TARGETS WITH MANY DIB DETECTIONS 



LSR velocity (km s ) 
-650 -600 -550 -500 -450 -400 -350 -300 -250 -200 -150 -100 -50 



Sightline 


3944.71 


3945.35 


4412.17 


MWISM 


* 


* 


Sp. type 


09.71b 


B6Ia 


B2.5 la 


avg. 


B0 la 


B9Ia 


Eb-v 


0.32 


0.33 


0.43 


0.32 








W(4428) 


1129 (106) 


628 (135) 


b 


865 


640 


370 


W(5705) 


< 31 


42 (25) 


< 19 


38 


5 


8 


W{5 /80) 


318 (25) 


202 (23) 


190 (1 /) 


163 


2 


30 


W(5797) 


62 (35) 


101 (28) 


45 (21) 


49 


4 


9 


W(6203) 


159 (65) 


93 (54) 


92 (40) 


54 


5 


4 


W(6269) 


60 (23) 


b 


< 14 







7 


W(6283) 


671 (75) 


587 (73) 


491 (82) 


362 


1 


10 


W(6379) 


b 


25 (18) 


63 (13) b 


28 


8 


1 


W(6613) 


80 (25) 


123 (20) 


57(15) 


68 





5 


W(6660) 


< 25 


54 (20) 


< 15 


16 





1 


W(6993) 


45 (23) 


77 (19) 


< 14 


38 





2 



Note. — DIB equivalent widths are in mA. Average DIB strengths (per 
unit Eb-v) for the Galactic diffuse ISM (scaled t o Eb-v = 0.32) are given in 
the column labelle d MW ISM (data from iHerb ig 1993. ThorburnJeTay|251 
Megier et al. 2005 and Cordinei 2006). No large sample of A6269 measurements 
were available. Possible stellar line contributions (in mA) to DIB equivalent 
widths are given in the columns labelled * for B0 la and B9 la stars, measured 
from synthetic stellar spectra with temperatures of 28,000 and 10,000 K, respec- 
tively. 

b Uncertain (or unmeasurable) equivalent widths due to line blends. 

DIB originated in the MW halo, and inconsistent with the 
values o bserved previously in the Milky Way ISM (see Fig- 
ure 8 of iMcCall et al.ll2010t) : 3) Assuming the Galactic halo 
obeys the average MW DIB W vs. E B - V relation, the equiv- 
alent widths expected from the M 3 1 foreground reddening 
of E B - V = 0.06 ± 0.02 are only 32 ± 11 mA for ,15780 and 
71 + 24 mA for /16283, whereas the measured DIB equivalent 
widths towards 4412.17 of 190 ± 17 mA and 491 ± 82 mA, 
respectively, are consistent with the M31 reddening (E™ :3 y) 
in this sightline; 4) The Nai spectrum of 4406.87 in Figure 
[3] which shows no evidence for M 3 1 Na i absorption, gives 
another indication of the properties of the foreground (MW) 
neutral gas in Field 3. The quantity and velocity of this fore- 
ground gas is very similar to the amount of gas in the other 
fields where the MW /15780 components (which can be easily 
separated from the M 3 1 DIBs), are clearly very weak. As- 
suming the foreground gas chemical properties are not grossly 
variable on the scale of the ~ 1 degree that separates the 
GMOS fields, this strongly suggests that the foreground DIBs 
in Field 3 would be similarly as weak. Small-scale structure 
in the Galactic halo could produce variations in the relative 
abundances of atoms, dust and DIB carriers, but given the 
reasons outlined above, it is highly unlikely that the Field 3 
DIBs originate in chemically peculiar, compact, high-velocity 
Galactic clouds. 

3.6. Hi 21 cm emission and the M 31 gas-to-dust ratio 

Dr. Robert Braun kindly made available data cubes 
and opacity - correc ted Hi column density maps from the 
iBraun et al.i j2009h 21 cm survey of M31, performed us- 
ing the Westerbork Synthesis Radio Telescope and the Green 
Bank Telescope. These data are at a spatial resolution of 18"x 
15", which corresponds to a linear resolution of 70 x 56 pc at 
the distance of M 3 1 . The H i spectra have a velocity resolu- 
tion of 2.3 kms -1 . From these data we used bi-linear inter- 







-650 -600 -550 -500 -450 -400 -350 -300 -250 -200 -150 -100 -50 
LSR velocity (km s" 1 ) 



Fig. 14. — Hi emission spectra derived from the survey by Braun et all 
120091) . M31 Nai velocities are plotted with dotted vertical lines. /U5780 
and 6283 DIB velocities (and \tr error bars) are plotted with black circles 
and red triangles, respectively. 
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polation to extract opacity-corrected column densities and H i 
spectra for our targets (shown in Figure [Pill. 

Opacity-corrected H i column densities and gas-to-dust ra- 
tios G = jiV(H ij/E^y are given in Table |4j the gas-to-dust 
ratio has thus been calculated assuming that half of the to- 
tal Hi column lies in front of the target stars. However, this 
fraction is uncertain, so the errors on G are ±100%. 

The average gas-to-dust ratios are 13.4, 8.2 and 7.7 (in units 
of 10 21 cm~ 2 mag -1 ), for Fields 1-3, respectively, with an 
overall average of 9.8. Previous measurements of the M31 
gas-to-dust ratio lie in the range «;5-15 x 10 2 1 cm 2 mag 1 
(Ivan den Berghl [T97l iBaiaia & Gergelvl[T977t iKumari [19791 
iLequeuxl 120001) . These values may be compa red with the 
nearb y Galactic value of 5.8 x 10 21 cirT 2 mag 1 dBohlin et al.l 
[19781) . which was derived from H i absorption lines and there- 
fore does not suffer from the problem of background contam- 
ination. It is difficult to draw firm conclusions from these data 
due to the uncertainties in the fraction of background H i in 
M 31, but is seems reasonable to conclude that the gas-to-dust 
ratio of M 3 1 is broadly similar to that of the Milky Way. Al- 
though many of the observed sightlines appear to have gas- 
to-dust ratios significantly greater than Galactic, this may be 
caused by an observational bias towards stars located on the 
near side of the M 3 1 disc; our derived gas-to-dust ratio would 
be an overestimate in such cases. Accounting for errors, the 
gas-to-dust ratio is lower towards 4412.17 than the average 
MW value. This sightline also clearly exhibits relatively weak 
DIBs per unit reddening. 

Fig.rBlshows that the velocities of the interstellar clouds as 
probed using Na i lines correspond imperfectly with the peak 
Hi velocity in many cases. This may be due to the presence 
of background H i clouds that do not contribute to the Na i ab- 
sorption spectra. However, for three of the Field 3 stars, the 
Na i velocities do not match the H i gas at all, with Na i ve- 
locities at least 20 kms -1 redwards of the Hi. For 4412.17 
and 4417.80, the Nai velocitie s are similarly displaced from 
the DIB velocities (see Section [3~5l ). but the DIBs match very 
closely with the positions of the Hi peak velo cities. Given 
that DIB carriers, Na i and H i are corr elated (e.g.lHobbsll974t 
iFerlet et alj|1985h lHerbidll993h iMcCall et alj|2010l) . and their 
velocities are generally closely matched in the diffuse ISM, 
there is likely to be a problem with the measured M 3 1 Na i 
cloud velocities in Field 3 (probably for t he reasons related to 
sky-subtraction discussed in Section [331 ). The measured Nai 
velocities would be offset to the red (as observed) if the MW 
component is not as strong as was derived in the vapid model 
fits or if the sky-line subtraction was incomplete. To resolve 
this situation would require high resolution Na i spectra in or- 
der to separate the M31 and MW components. The Hi gas 
in Field 3 spans a relatively narrow range of radial velocities 
(from about -120 to -60 kms~'), and should be considered 
as a more reliable indicator of the velocity of the interstellar 
matter in these sightlines. The resulting velocities from Gaus- 
sian fits to the H i emission profiles are shown for comparison 
with the M 31 DIBs in Field 3 in Figures[3to[II] 

CO peak velocities from radio observations bv lNieten et alj 
( 120061) are consistent with those of the Hi gas. For Field 3, 
CO matches better with H i than Na i, which provides further 
evidence that our measured M 31 Nai velocities for 4408.04, 
4412.17 and 4417.80 are probably in error. 

3.7. M31 H\\ regions 

In addition to the stars in TableQ] three H n regions were ob- 
served serendipitously in Field 3, spatially coincident with the 
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Fig. 15. — Sky-subtracted spectra from observed Hn emission regions in 
the telluric rest frame. The Nn 6583 A and H a lines used to derive oxygen 
abundances are marked. 

GMOS slits of the stars 4408.79, 4408.94 and 4417. 12. Using 
the relative strengths of the emission lines in these spectra, it 
is possible to derive an estimate for the interstellar metallicity 
in the vicinity of Field 3. The usefulness of these spectra is 
limited because they are not flux calibrated and the blue and 
red wavelength regions were recorded in separate exposures. 
Nevertheless, an estimate of the O/H abundance ratio can be 
obtained from the rel ative N n 6583 A and Ha line fluxes (see 
IPettini & Pajell2004h . 

The emission spectra were extracted from regions at the 
edges of the GMOS slits where the stellar flux was negligible, 
sky-subtracted using sky spectra from nearby slits showing no 
nebular emission, and continuum-subtracted using low-order 
polynomial fits. The final spectra are shown in Figure[l5] The 
logarithmic flux ratios (log([Nn]/Ha)) for 4408.79, 4408.94 
and 4417.12 are -0.54, -0.63 and -0.61 respectively, which 
correspond to oxygen abundances 12 + log[0/H] = 8.59, 8.53 
an d 8.51 (all with lcr er rors of +0.18), using equation (2) 
of IPettini & Pagell (120041) . The average value, 8.54 + 0.18, 
is similar to that in the solar neighbourhood, where the oxy- 
gen abundance of H n regions is in the range 8.45 - 8.50 (see 
iPilyugin et al.ll2003l and references therein). 

Our derived [O/H] ratios add to the limited data al- 
ready availab l e in th e literature on the metallicity of M 3 1 . 
iTrundle et alj d2002l) reviewed previous work on oxygen 
abundance measurements of H n regions in M 3 1 and iden- 
tified that there is some uncertainty in the metallicity gra- 
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dient across M31. Measured oxygen abundance gradients 
range from -0.013 to -0.027 dexkpc -1 , depending on the 
method used to derive [O/H]. Extrapolating these gradients 
gives 12+log[0/H] = 8.70 to 9.21 at the centre of M 3 1 , com- 
pared with 8.90 in the centre of the Milky Way. The three H n 
regions we observed are approximately 23' from the centre of 
M31, which corr esponds to 5.2 kpc (assuming a distance of 
783 kpctoM31;|Holland|[l99!. The average oxygen abun- 
dance in the three Hn regions (8.54), lies towards the lower 
end of the range of values previously reported at this radius in 
M31. 

Fields 1 and 2 lie at radial distances (r) of 63' (14.3 kpc) 
and 41' (9.3 kpc) respectively from the centre of M31. Us- 
ing the metallicity relation of 12 + log[0/H] = 8.7 - 0.013r, 
which was deduced from the study of B star abundances by 
iTrundleet aD (120021) . we find oxygen abundances of 8.51 and 
8.58 for these two fields. Thus, there is no evidence to sug- 
gest that the metallicity of the ISM along any of our surveyed 
sightlines is significantly different from solar. 

3.8. Interstellar radiation field 

Near-ultraviolet interstellar radiation field strengths (NUV 
ISRFs) were calculated for the M 3 1 ISM in each sightline, 
using the NUV GALEX mosaic of M 3 1 (imaged at a cen- 
tral wavelength of 2271 A) from iThilker et al.1 (I2005 S). The 
NUV ISRF strengths (Inu v) are given in Table |4] in units of 
the Galactic NUV ISRF of lDrainel(ll978l) . integrated over the 
aperture size and NUV wavelength response function of the 
GALEX telescope. Maps of Inuv for each of our three fields 
were computed from the GALEX image and are shown in Fig- 
ure E] 

To calculate these Inuv maps, individual pixel fluxes were 
background-subtracted and corre cted for foreground ex tinc- 
tion using Anuv = 8.0 x Eb-v dGil de Paz et al.l 120071) . as- 
suming a MW foreground extinction of 0.06 mag (see Sec- 
tion [3T4ji. The M31 extinction was derived for each field 
from the average of the E^v values of the stars listed in Ta- 
ble Q] For each sightline in the ISRF maps, Inuv was calcu- 
lated at a point at a line-of-sight distance r halfway into the 
foreground M 31 cloud, using th e extinction-distance relation 
derived by IVergely et al.l (119981) for the solar neighbourhood 
(dA v /dr = 1.2 magkpc -1 , with R v = A v /E B -v = 3.1), and 
assuming that the observed NUV flux is emitted from the im- 
age plane (at a radial distance of 783 kpc). The NUV photons 
were assumed to be subject to the same level of extinction as 
a function of distance as in the solar neighbourhood; conse- 
quently, the calculated Inuv are dominated by the flux from 
OB stars in the immediate vicinity of the interstellar clouds 
of interest; stars at distances greater than 1 kpc make a neg- 
ligible contribution to the total NUV flux at a given point in 
space. The NUV ISRF strengths thus calculated are only ap- 
proximate because they take no account of the fact that the 
shape of the extinction law, and the extinction as a function of 
distance in M 31, may be different from that observed in the 
Milky Way. The M 3 1 foreground extinction is treated as a 
uniform slab of material, so any small-scale inhomogeneities 
in the M 3 1 ISM are neglected. 

Within Fields 1-3 respectively, the interstellar UV field 
strengths (Inuv) averaged over all pixels are 0.24, 0.77 and 
0.34 (in units of the Draine field), with an average value of 
0.45. Evidently, even though these fields are relatively rich in 
OB stars, the NUV ISRF is still significantly weaker than that 
of the Milky Way. This is consistent with the conclusions of 



previous authors (see Section[TJ. Field 2 has the greatest Inuv 
values due to the presence of OB78 (NGC206) - the richest 
OB association in M31. Indeed, the sightlines towards the 
DIB target stars in this field have NU V field strengt hs which 
are comparable to or greater than the dDraindl 19781) Galactic 
value. On the other hand, most of the stars in Field 1 to- 
wards which DIBs are observed have Inuv < 0.5, including 
the star in our sample 3944.71, towards which the M31 DIBs 
are strongest. 

3.9. Aromatics and dust emission 

For the three M31 fields presented here, we obtained 
reduced infrared Spitzer Space Telescope images in the 
IRAC 8 urn and M IPS 24 urn bands from the surveys of 
iBarmbv et all ( 120061) and lGordon et al.1 ( 120061) . The flux in the 
8 /mi band is dominated by PAH emission fro m the 7.7 /an 
(and 8.6 /an) aromatic band s (see, for example iDraine &~Q 
1200 11 and iFlagey et aPl2006l) . whereas the flux in the 24 /an 
band predominantly arises from interstellar dust. The ratio of 
8 /mi-to-24 /mi fluxes may therefore be used as an approxi- 
mate measure o f the relative PAH emission intensity per unit 
dust mass (e g. IFlagey et al.1 l2006t lEngelbracht et al.l 120081: 
iGordon et al.l 120081) . and as such, provides an indication of 
the abundance and excitation level of PAHs in the ISM. 

The 8 //m IRAC images are shown for our three fields 
in Figure [17] The interstellar emission flux was measured 
for each of the sightlines in Table [4] from the background- 
subtracted Spitzer images using a ring median filter of radius 
5 pixels and width 2 pixels in order to exclude the flux from 
stellar point-sources. To measure the 8/24 //m interstellar flux 
ratio, the 8 /an images were convolved and rebinned to the 
lower spatial resolution o f the 24 /an images u sing the convo- 
lution kernel provided bv lGordon et al.l (120081) . 

In the plane of the Milky Way, the interstellar 8/24 /mi flux 
ratio ( R) varies betw een about 0.1 and 2.0, with a mean value 
of 1.6 (Fl ageyl 120071) . Similarly, in nearby spiral galaxies, R 
also varies between about 0. 1 and 2, with consid erable scat- 
ter within individual galaxies (B endo et al.ll2.Q08b . The ratio 
R correlates with metallicity, which has been interpreted as 
being due to variation in the strength of PAH emission rel- 
ative to the 24 /an dust emission, at least partly as a resul t 
of differences in PAH abundances (lEngelbracht et al.l 120051) . 
However, PAH emission co rrelates better with i onization of 
the IS M than the metallicity (IGordon et al.ll2008l;lBendo et al.l 
2008), so the 8 /mi flux may be a better indicator of the degree 
to which the PAHs are ionised rather than their overall abun- 
dance. R also tends to be lower inside H n regions and higher 
in denser, dustier gas associated with PDR boundaries. 

In our GMOS fields, R values range between 0.63 and 1.78, 
with mean values of 1.05, 1.16 and 1.04 in Fields 1-3, respec- 
tively. Thus, there is nothing particularly unusual about the 
PAH emission in these fields compared with other galaxies, 
although the mean values are somewhat lower than observed 
in the Milky Way. This is perhaps indicative of a reduced 
level of ionisation or processing of the aromatic material in 
our M 3 1 fields as a result of the lower UV radiation fields 
that are typical of M 3 1 . 

The strength of the 8 /mi emission discussed here is at vari- 
ance with previous studies (Cesarsky et al. 1998; Pagani et al. 
119991) that identified anomalously weak aromatic mid-IR 
emission bands in the bulge and nucleus of M 3 1 . Thus, the 
PAH emission traced by the 8 /mi flux, at least in our observed 
fields (which are away from the nucleus/bulge), is apparently 
not as weak as previously thought, and implies that the PAHs 
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in M 3 1 may in fact not be dissimilar to those found in the 
Milky Way and other spiral galaxies. It must be noted, how- 
ever, that the NUV ISRF in our observed fields may be greater 
than is typical for M 3 1 due to the concentration of OB stars, 
which could give rise to increased 8 fim emission in these re- 
gions relative to the average conditions in M 3 1 . 

4. DIB STRENGTHS AND CORRELATIONS 

Figure [16] shows the relationships between the equivalent 
widths of a selection of the strongest diffuse interstellar bands 
(X14428, 5780, 6283, 6379, 6613 and 6660) and interstellar 
reddening CE^_y), for sightlines in M 31 and the Milky Way. 
All of the newly-observed M 3 1 DIB measurements lie within 
the range of values previously observed in the Milky Way, 
apart from /16660, which seems anomalously strong, perhaps 
due to spectral contamination, which may be indicated by ref- 
erence to the telluric standard spectrum in Figure fT2l 

The majority of the M 3 1 data points plotted in Figure [16] 
lie within the middle-to-upper part of the range of values cov- 
ered by the MW data. Thus, per unit E B -v, the M 31 DIBs 
are as strong or stronger than those typically observed in the 
Milky Way. The relative enhancement of DIB strengths is 
most pronounced for Field 1 (shown as black filled circles in 
the plots), for which the equivalent widths of the strongest 
(and therefore most reliably measured) DIBs /M5780, 6283 
and 6613 are clearly at the upper end of the range of Galac- 
tic values. The locations of those sightlines with the strongest 
DIBs are shown in Figure [17] the sizes of the circles indi- 
cate the relative strength of the /16283 DIB per unit Eb-v in 
each sightline. It is evident that the ISM in Field 1 hosts the 
strongest DIBs per unit Eb-v and Fields 2 and 3 the weakest. 
The DIBs with the largest equivalent widths in our sample are 
towards 3944.71, 3945.35 in Field 1. These are the sightlines 
with the most reddening and the greatest Hi column densi- 
ties. Relative to the amount of dust in the 3944.71 sightline, 
the /1/15780, 6283 and 6613 equivalent widths are at the very 
upper limit of what has previously been observed in the Milky 
Way, so these sightlines have exceptionally strong DIBs. As 
can be seen in Table[3] the majority of the measured DIBs in 
the sightlines towards 3944.71 and 3945.35 are significantly 
stronger than those in the Milky Way, for the same reddening. 
These two sightlines have relatively low NUV ISRF strengths 
- about 0.5 and 0.4 times th at of the Milky Way, respectively. 

As derived in Section 13.41 the mean /15780 equiva- 
lent width per unit reddening measured in the Milky Way 
is [W(51S.0)/E B - V ] mean = 0.53 + 0.23 Amag" 1 . Using 
data on the /16283 DIB from the same Galactic literature, 
[W (622,3) /Es-vlnean = 1-18 + 0.53 Amag- 1 . The errors are 
the standard deviations of the respective data. The mean val- 
ues for W(5780)/£Vv in Fields 1, 2 and 3 are 0.93, 0.7 and 
0.73 Amag- 1 , respectively, and for W(62S3)/E B -v, 2.4, 1.37 
and 1.27. These values are consistent with the identified trend 
that our observed M31 DIBs are relatively strong compared 
with the Milky Way, and highlight in particular, the strengths 
of these DIBs in Field 1. 

The weakest DIBs (per unit Eb-v) are observed towards 
4030.28, 4034.61 and 4037.92 in Field 2 and 4412.17 in 
Field 3. These four sightlines are near to regions of strong 
UV flux (see Figure ITTb. and include the strongest NUV IS- 
RFs of the observed sightlines, with Inuv Z 1- The small- 
est W(6283)/£V\/ value is observed towards 4030.28, where 
Inuv is largest, and the largest W(6283)/£Vy value is ob- 
served towards 3936.51, where Inuv is smallest, which is 
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Fig. 16.— AA442&, 5780, 6283, 6379, 6613 and 6660 DIB equivalent 
widths as a function of Eb-v for M31 and the Milky Way. Data points 
for t he three GMOS fields are plotted in different colours. Measurements 
from [Cordiner et al. 1 2008a) are included (blue squares). DIB equivalent 
width upper limits are shown with downward-po i nting arrows. MW DIB 
data are f r omlHerbid rf99lh . iTho"rburn et alj 1 2003), Megier etajj <2005l) and 
ICordineJ 120061) . and for i4430. IHerbieH 19751) . For the AA5180 and 6283 
DIBs, linear regressions are shown for the MW data (black solid lines) and 
the M 31 data (blue solid lines). A linear fit that takes into account the errors 
in both ordinates is also shown for the M31 /1/15780 and 6283 data (dashed 
blue lines). 
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Fig. 17. — Three-c olour composite maps of our observed M31 fields. Blue corresponds to th e NUV ISRF (see Sect. 13. 8t . green is the optical V-band image 
i Massev et al. 2006), and red is the Spitzer IRAC 8 fim data that traces PAH emission (Sect. l3~9l . Intensities are shown on a logarithmic scale. ISRF values range 
from 0.05 to 0.5, 6.0 and 1 .0 Draine for Fields 1-3, respectively. 8 /im flux values range from 0.0 to 50.0 MJy sr~' for all fields. Lines-of-sight with M 3 1 DIBs 
are marked with circles, the sizes of which denote the strength of the /16283 DIB per unit reddening (relative to the mean M31 DIB strength): small, 10"-circles 
denote weak; medium, 20"-circles denote average and large, 30"-circles denote strong DIBs. Crosses indicate sightlines with no DIB detections, squares denote 
upper limits. Also marked are the tell uric standards (t ) , the H n emission regions (e) and the single Wolf-Rayet star observed in Field 1 (wr). The two targets with 
previously observed DIBs reported by Cordiner et al. (2008a) are about 6 arcminutes (s» 1.3 kpc) south of Field 2. 



suggestive of a correlation between the interstellar UV field 
and DIB strengths. To analyse this correlation, the equivalent 
widths per unit reddening of the two strongest DIBs (with the 
smallest relative measurement errors), /L15780 and 6283, are 
plotted with respect to the NUV ISRF strength in Figure [T8l 
There is a significant negative correlation between W/E^ v 
and Inuv f° r both DIBs at the 68% confidence level, shown 
by the grey polygon. The correlation coefficients between the 
data for /U5780 and 6283 are -0.55 and -0.41, respectively. 
This indicates that as the interstellar UV flux increases, the 
abundances of both of these DIB carriers tend to become less, 
relative to the amount of dust. Plausible explanations for this 
phenomenon are that the DIB carriers become ionised or are 
photodissociated, or that the chemical environments in which 
the carriers form are destroyed by NUV radiation. 

A link between UV flux and DIB strengths is corroborated 
by the pattern of DIB observations in Field 2 and the mor- 
phology of the ISM in this region. This field covers the OB78 
association, which is the largest, most active region of recent 



star formation in M 3 1 . Most of the stars at the heart of this 
association do not show significant reddening or DIBs (see 
Figure [17] ) - the majority of the gas and dust around these 
stars has been blown away by strong stellar winds and su- 
pernovae dBrinks & Shandll984t lLoinard et alj|1996l) . creat- 
ing a large, hollow region in the ISM several hundred parsecs 
across. There is still sufficient interstellar matter at the periph- 
ery of the central cluster that DIBs are observed towards many 
of the nearby stars, and the 8 fim flux is intensely strong as a 
result of the high UV-optical flux that excites the interstellar 
aromatics. The DIBs are generally quite weak in this field 
compared to the M 3 1 average, which we suggest is due to the 
destruction by strong UV irradiation of the environments in 
which the carriers form. 

The relationship between DIB strengths and UV flux has 
been examined previously in the Milky Way and in the Mag- 
ellanic Clouds. Indeed, the relative weakness of the DIBs ob- 
served in the Magellanic Clouds may be attribut able to the 
strong radiation fields that pervade these galaxies (ICox et al] 
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2006: C ordinerj l2006t iWeltv et al.l 120061) . UV-rich regions 
such as those found in the vicinity of the Orion Nebula and 
Scorp ius are also known to show weak DIBs dCami et al.l 
U997tlSnowet al.lll995l) . IWeltv et all (120061) found no signif- 
icant correlation between UV flux and DIB strengths in the 
MCs. However, the UV field strengths in their analysis were 
derived from H2 excitation, and it is known that the DIB carri- 
ers do not predominantly resi de in the relati vely dense regions 
of the ISM traced by H 2 (e.g. lHerbig|[T993l) . 

The narrow (AA5797, 6613, 6660 and 6993) DIBs are en- 
hanced towards 3945.35 with respect to the average MW 
diffuse ISM, while the broad DIBs are not. On the other 
hand, the line of sight towards 3944.71 (where the DIBs are 
strongest), shows enhanced broad (/i/15780, 6203 and 6283) 
DIBs with respect to the a verage MW value. Previous stud- 
ies in the Milky Way (e.g. ICami et al.ll!997t iKrelowski et alJ 
119981) . identified that these broad DIBs tend to be strong in 
relatively strongly UV-irradiated sightlines (where the narrow 
DIBs are weak), whereas narrow DIBs are strong in well- 
shielded sightlines. However, the negative correlation we 
have identified between the strengths of the /L15780 and 6283 
DIBs and Inuv suggests that the environments or the carri- 
ers of these DIBs are destroyed by UV radiation. Their large 
strengths in relatively strongly-irradiated Galactic sightlines 
seems to be at odds with our observed correlation, but could 
instead be evidence for the production of enhanced quanti- 
ties of DIB carriers as a result of photochemistry or dust 
photo-abrasion in these Galactic regions. Indeed, it is notable 
from Figure[T7]that the two sightlines with the strongest DIBs 
(3944.71 and 3945.35), reside in a region of both strong 8 fim 
flux (indicative of dense, dusty gas) and moderately high UV 
field strength. Given the ass ociation of mid-IR emission with 
ionised PAHs dTielensl2008l) . this suggests that the carriers of 
the /W5780, 6203 and 6283 diffuse bands may well be ionised 
molecules. 

The correlations between /M5780 and 6283 and Eb-v are 
shown in Figure Q~6] Taking into account the error bars on 
W and Eb-v, least-squares linear fits to the data result in the 
parameters shown in Table [6] Within the quoted l<x errors, 
there are no significant differences between the slopes and in- 
tercepts of M31 and MW data for either DIB. Kolmogorov- 
Smirnov (K-S) tests were applied to the residuals of the ob- 
served data compared with the MW fits. The resulting proba- 
bilities (1 - P) describe the likelihood that the M 31 and MW 
DIB data have different distributions with respect to Eb-v- For 
/15780, (1 — P) = 0.93, which shows there is a high statisti- 
cal probability that the M 31 and MW DIBs follow a differ- 
ent relationship with Eb-v- This is consistent with our ob- 
servation that the M 3 1 DIBs are typically rather strong com- 
pared with in the Milky Way. For /i6283, the derived value 
of (1 — P) = 0.50 is inconclusive, and indicates that more 
observational data are needed. 

Relatively weak DIBs are observed towards 4412.17 which 
has a very low G value; relatively strong DIBs are observed 
towards 3936.51, which has a high G value. These results 
could be explained if variations in G are caused by the de- 
struction of dust grains (by shocks and UV radiation), and 
if the grain degradation products contain or give rise to the 
DIB carriers. Conversely, if gas-phase species are depleted 
out onto the dust, reducing the gas-to-dust ratio, then molec- 
ular DIB carriers would also accrete onto the grains, making 
the bands weaker. However, when accounting for the sub- 
stantial errors on G across all sightlines, there is no statisti- 
cally significant correlation between the gas-to-dust ratio and 



TABLE 6 

Fit parameters for M 31 /l/i5780 and 6283 DIB equivalent widths 

VS. REDDEMNG 



DIB intercept a slope b correlation coefficient R 

MW5780 14 ± 10 469 ±29 0.81 

M 31 5780 73 ±62 426 ± 234 0.60 

MW6283 123 ± 80 876 ± 174 0.55 

M31 6283 146 ± 102 1000 ± 430 0.57 



Note. — Equiv alent width W = a + (b X Eb-v)- The outlier data 
from MAG 70817 (Cordiner et al. 2 0083) were omitted from the fits. 



W(5780) or W(6283), as shown by the grey polygons in Fig- 
ure [18] We know of no prior reports of correlations between 
DIB strengths and G in the MW; in the LMC, W/E B -v and G 
have been shown to be uncorrected, albeit with a very limited 
number of data points (ICox et al.l2006l;ICordinerl2006l) . Thus, 
an alternative explanation may be required for the weakness 
of the DIBs towards 4412.17, perhaps related to the relatively 
high strength of the 8 fim emission at the position of this 
star. This indicates the presence of dense, dusty gas of the 
kind where DIBs t end to be relatively weak in the Milky Way 
(lSnowetal.1119951) . 

In Section r3~9l the relationship between the mid-IR 8/24 mi- 
cron Spitzer flux ratios and the abundance of PAHs was dis- 
cussed. Given that some of the DIBs are hypothesised to be 
caused by electronic transitions of neutral and ionised PAHs 
(the same classes of molecules believed to cause the 8 fim 
emission), it is of interest to look for a correlation between 
the DIB strengths and the mid-IR flux ratio. However, as 
shown by Figure [18] there is no relation evident between the 
strengths of /M5780 & 6283 per unit Eb-v and the 8/24 mi- 
cron ratios observed in our M 3 1 sightlines. This indicates that 
the carriers of these particular DIBs are not the same as the 
aromatics which are predominantly responsible for the 8 /zm 
emission. There is some uncertainty in this result due to the 
fact that - in common with the Hi 21 cm emission - it is im- 
possible to separate the background contributions to the IR 
flux in our observed sightlines, so no definitive conclusions 
may be drawn. Nevertheless, it is apparent from the maps in 
Figure \TT} that the distributions of sightlines with strong or 
weak DIBs do not generally follow the 8 fim emission inten- 
sity. 

The exc eptionally stro n g DIBs that were previously de- 
tected by iCordiner etail d2008al) towards MAG 70817 in 
OB78 now seem to be unusual compared with typical M31 
and MW trends. The strong /L5780 DIBs rep orted in the more 
extensive Keck DEIMOS study of OB78 bv lCox & Cordiner! 

were known to have large error bars on W and Eb-v, 
and the deduction of that study, that the M 3 1 DIBs are sub- 
stantially stronger than those observed in the Milky Way, is 
not supported by the more accurate and statistically meaning- 
ful results presented here. 

Past studies, highlighted in Section [T] have discovered dif- 
ferent dust extinction and polarisation properties in M31 
compared to the MW, and in particular, a lack of small 
graphitic dust grains and an abundance of (relatively less UV- 
processed) PAHs. The similarity of the DIB spectrum to 
that in the MW indicates that the cause of these differences 
in the dust properties has little impact on the DIB carriers. 
That markedly different UV-extinction and optical polarisa- 
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Fig. 18. — Correlations between ,1/15780 and 6283 DIB equivalent widths per unit Eg-v and interstellar physical / chemical properties: Interstellar UV field 
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excludes blended DIBs and upper limits. 



tion curves have been observed in M 3 1 adds to the body of 
evidence that DIB carriers are not closely related to larger in- 
terstellar dust particles or to the smaller graphitic grains be- 
lieved to be responsible for the 2175 A UV extinction bump. 
However, studies of M 3 1 dust and extinction have so far been 
very limited, and little is presently known about the detailed 
characteristics of the average M31 extinction curve. Further 
dedicated studies of this are warranted. 

The general lack of UV/FUV radiation throughout M31 
compared with the Milky Way (see Section [TJ, due to the rel- 
atively low numbers of young stars, apparently does not have 
a detrimental effect on the DIB carrier abundances. Thus it 
seems that the processing of interstellar material by strong, 
hard UV radiation may not be required to produce DIB carri- 
ers. There is still, however, a significant NUV flux throughout 
much of M 3 1 (which, in the vicinity of dense clusters of OB 
stars, matches or exceeds the MW average), so the require- 
ment for lower-energy UV photons in the production of DIB 
carriers cannot be ruled out. 

5. CONCLUSION 

We obtained optical spectra of over 30 early-type stars in 
three fields in M31 using the Gemini Multi-Object Spec- 
trograph. Spectral classifications were presented, compris- 
ing the first spectral types for 20 of these stars. Redden- 
ings (Eb-v) towards our targets were calculated, including the 
foreground-corrected M 3 1 component. The M 3 1 and MW 



components of interstellar Nai towards the most-reddened 
subset of our targets were measured and their their column 
densities, Doppler widths and peak radial velocities derived. 

Eleven diffuse interstellar bands (/U4428, 5705, 5780, 
5797, 6203, 6269, 6283, 6379, 6613, 6660, and 6993) have 
been detected in the M 3 1 ISM, with profiles and patterns of 
strengths similar to those found in the Milky Way and with 
radial velocities generally matching the interstellar Nai and 
Hi and the stellar velocities. Measured radial velocities of 
the stars and interstellar matter match those of previously ob- 
served emission-line objects at the same projected galactocen- 
tric radii in M 3 1 . 

The M 31 /U5780 and 6283 DIB equivalent widths as a 
function of Eb-v fall within the envelope of values typically 
seen in the Milky Way. However, contrary to LMC and SMC 
DIBs which are weaker than those in the Galaxy, the M31 
sightlines generally have DIBs that are as strong or slightly 
stronger than the average Galactic values. The M31 and 
MW DIB W vs. Eb-v distributions were compared using a 
K-S test, which yielded a 93% probability that the observed 
M31 /15780 DIB data are different from that of the MW. 
However, this statistic does not account for the errors in the 
data values. We measure a mean M 31 /15780 DIB equivalent 
width per unit reddening of [W(51S0)/E B - V ] mea „ = 0.78^, 

compared with 0.53 + 0.23 Amag~' in the MW. Similarly, 
[W (6283) /E B -v]m e an = 1.81 i^Qgj A mag -1 , compared with the 
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MW value 1.18 ± 0.53 A mag 1 . Thus, although the ob- 
served mean values are indicative of stronger-than-average 
DIBs within M31, the sizes of the error bars preclude the 
identification of any significant differences between the DIB 
strengths in M 3 1 and the Milky Way. Higher signal-to-noise 
observations of more M 3 1 sightlines are required to resolve 
this issue. 

A correlation between DIB strengths and Inuv has been 
found at a 68% confidence level, which indicates that the 
DIB carrier abundances are enhanced in regions with weaker 
UV interstellar photon flux. We therefore hypothesise that 
the general lack of U V photons in parts of M 3 1 benefits the 
survival of DIB carriers and their chemical precursors against 
photo-destruction, and that FUV radiation is not a require- 
ment for the production of DIB carriers. 

The DIBs observed in the spectrum of 3944.71 are very 
strong per unit Eb-v- It is clear, however, tha t the unusually 
strong DIBs previously detected in M31 by iCordiner et al.l 
(2008a]) and ICox & Cordinerl d2008b are anomalous with re- 
spect to the average M 3 1 trend, and that in general, the rela- 
tionship between M 3 1 DIB equivalent widths and reddening 
is very similar to that of the Milky Way. The degree of scat- 
ter of the M 3 1 DIBs around the average relationship (per unit 
Eb-v), is also similar to that observed in the Galaxy. 

The similarity of the M 31 DIB W/E B - V relation to that of 
the Milky Way suggests that DIB equivalent widths may be 
applicable as probes of Eb-v in other, similar spiral galax- 
ies. However, there is evidence that the relationship for 
M 3 1 is slightly different, perhaps as a result of the low UV 
field strength, and therefore further study of DIBs in exter- 
nal galaxies is warranted. Particular caution must be exer- 
cised when using DIBs as a measure of reddening in galax- 
ies with substantially different metallicities, gas-to-dust ratios 
and interstellar radiation field strengths (such as the Magel- 
lanic Clouds). 

When taking into account the substantial errors on G, 
the strengths of the /U5780 and 6283 DIBs (measured by 
W/Eb-v) do not show any significant correlation with the gas- 
to-dust ratio G. 
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No significant correlation is found between W/Eb-v and the 
8/24 /im flux ratio for the AA57W and 6283 DIBs, indicat- 
ing that the (PAH) carriers of the 8 fim infrared band are not 
closely related to these DIB carriers. The 8/24 fim and gas-to- 
dust ratios in our fields are in the range of values previously 
observed in the Milky Way, so in these respects the gas and 
dust properties of M 3 1 are not found to be distinguishably 
different from those of the Galaxy. 

Finally, it should be noted that our results reflect the be- 
haviour of the ISM in only three small regions of M 3 1 and 
therefore do not necessarily reflect the global behaviour of 
the M 3 1 disc. However, our regions do span a range of galac- 
tocentric distances and have a range of different properties in 
terms of the structure and composition of the gas, the density 
and type of stellar objects, the radiation field strength and the 
metallicity, and are therefore indicative of the DIB behaviour 
expected throughout much of the rest of the disc. 
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